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Abstract. A matrix model of an asymptotically free theory with a bound state 
is solved using a perturbative similarity renormalization group for hamiltonians. 
An effective hamiltonian with a small width, calculated including the first three 
terms in the perturbative expansion, is projected on a small set of effective basis 
states. The resulting small hamiltonian matrix is diagonalized and the exact 
bound state energy is obtained with accuracy of order 10%. Then, a brief 
description and an elementary illustration are given for a related light-front 
Fock space operator method which aims at carrying out analogous steps for 
hamiltonians of QCD and other theories. 


1. INTRODUCTION 

This lecture has two aims. The first aim is to show a simple example of a 
new kind of calculation of effective hamiltonians, based on the perturbative 
similarity renormalization group. 0 i The second aim is to show how one can 
generalize the simple example and start systematic perturbative calculations 
for quantum field theoretic hamiltonians in the light-front Fock space. 

Although the methods we present are quite general, the main motivation 
came from QCD. QCD is asymptotically free and its perturbative running 
coupling constant grows at small momentum transfers beyond limits. This 
rise invalidates usual perturbative expansions in the region of scales where the 
bound states are formed. 

Ref.|^ outlined a light-front hamiltonian approach to this problem in QCD, 
using the perturbative similarity renormalization group. Independently, Weg¬ 
ner proposed a flow equation for hamiltonians in solid state physics. He 
introduced an explicit expression for the generator of the similarity transfor¬ 
mation which leads to a Caussian similarity factor of a uniform width. 

Wilson and I have solved numerically a simple matrix model to gain quanti¬ 
tative experience with the similarity scheme using Wegner’s equation. We 
also made perturbative studies. 0 This lecture is based on those works in the 
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part describing the model. The remaining part contains an outline of how one 
can attempt to make similar steps for light-front hamiltonians in quantum field 
theory using creation and annihilation operators. 


2. MODEL 

Consider a quantum theory which is characterized by a large range of energy 
scales as measured by certain Hq. QCD has this feature. It extends in energies 
from oo (asymptotic freedom) down to the infrared energy region. We represent 
the theory by a model with a hamiltonian H = Hq + Hi acting in a space 
spanned by a finite discrete set of nondegenerate eigenstates of the hamiltonian 
Ho, 

Ho\i > = E,\i > . (2.1) 

Matrix elements of the interaction are assumed to be 

<i\Hi\j> = -g^/E^Ej . (2.2) 

g is a dimensionless coupling constant. 

We choose Ei = 2* and M < i < N. M is large and negative and N is large 
and positive. We use M = —21 and = 16 in our numerical example. Let the 
energy equal 1 correspond to 1 GeV. Then, the ultraviolet cutoff corresponds 
to 65 TeV and the infrared cutoff corresponds to 0.5 eV. 

The same model can be alternatively derived by discretization of the 2- 
dimensional Schrodinger equation with a potential of the form a coupling con¬ 
stant times a (5-function. [|| 

For g > 1/38, the hamiltonian matrix has one negative eigenvalue and 37 
positive eigenvalues, g is adjusted to obtain the negative eigenvalue equal — 1 
GeV; g ^ 0.06. This eigenvalue corresponds to the s-wave bound state energy 
in the 2-dimensional Schrodinger equation. 

We calculate effective hamiltonians, Ti = using the similarity renor¬ 

malization group equations in the differential form. The effective hamiltonians 
are parametrized by their energy width A. The notion of the hamiltonian width 
will become clear shortly. We use Wegner’s flow equation Q 

= -y.m'nin], (2.5) 

with the initial condition Ti.{oo) = H. The matrix H is the diagonal part of Ti. 
with elements = TimmSmn- Thus, 7i(A) is a unitary transform of H and 
both have the same spectrum (see Wegner’s lecture in this volume). 

Equation (2.5) can be approximately solved for a small g keeping only terms 
order 1 and g. One obtains 

Timn = Em^mn “ 5\/ EmE^ exp [-{Em “ E^f /\^] . (2.6) 




RENORMALIZATION OF HAMILTONIANS 


3 


Here, Prnm = (1 ~ g)Em- The Gaussian factor of width A is the similar¬ 
ity function. This explains the notion of the hamiltonian width. Ref. 
demonstrated that the Wegner flow equation has a renormalization group in¬ 
terpretation. Including terms order g^, we let g depend on A and we introduce 
g{X) = g. It follows from equations satisfied by the matrix elements Timn with 
the indices m and n close to M that, neglecting small energies, 

dg/dX = -g^^'^exp[-2E^/X^] , (2.7) 

£ 

and 5 ( 00 ) = g. Analytic integration of Eq. (2.7) in the model gives, approxi¬ 
mately, 

ga{X) = (1.451ogA-0.9)-^ (2.8) 

5 a (A) grows when A gets smaller and it exhibits the asymptotic freedom behav¬ 
ior: it is smaller for more violent interactions (i.e. of wider range in energy). 
5 a (A) blows up to infinity for A ~ 1.9 GeV. In this approximation, matrix 
elements of H for Em ^ En ^ A can be written as 

'XCmniX') = Em^mn 9a{X) \/ EmE^ exp ^ [Em Eji\ /X ] -f COTTCCtioTlS . 

(2.9) 



Figure 1. The approximate running coupling 5 a (A) from Eq. (2.8) and the exact 
running coupling 5 (A), plotted as functions of the effective hamiltonian width A. The 
matrix element jl{X) = 7i-i,-i(A) — 0.5 GeV is also plotted to show the width range 
where the bound state eigenvalue appears on the diagonal. 
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Now, = [1 — ga{X)]Em. The energy order of low energy states is 

reversed when ga(A) grows above 1. 

The exact running coupling, g{X), is defined by writing Hm,m+i{X) = —g{X) 
^^EmEm+i- Eq. (2.9) shows that ga = g for large A. To find g for all values 
of A, we solved Eq. (2.5) numerically. Fig. 1. shows that the approximate 
solution blows up in the flow before the effective hamiltonian width is reduced 
to the scale where the bound state is formed. That scale, order 1 GeV, equals 
A at which the bound state eigenvalue appears on the diagonal. The diagonal 
matrix element is also shown in Fig. 1. 

The key feature, visible in Fig. 1, is that the exact effective coupling con¬ 
stant does not grow unlimitedly. The similarity renormalization group for 
hamiltonians provides a new option for investigating bound state dynamics in 
asymptotically free theories. The question is how far down in A we can reach 
using perturbation theory instead of the exact solution. The answer is: down 
to 1 GeV in second order with 10% accuracy. This is illustrated in Fig. 2. 

The remaining question of how small the space of states can be on which one 



Figure 2. The accuracy of the bound state eigenvalues obtained from effective 
hamiltonians whose renormalization group flow with the width A is calculated ex¬ 
panding in powers of the effective coupling constant g{Xo) and including terms order 
1, g(Xo) and g^{Xo). The accuracy is given as ratio of the bound state eigenvalue ob¬ 
tained by diagonalization of the effective hamiltonian of width A to the exact value, 
— 1 GeV. The curves correspond to the indicated values of Ao (in units of GeV). The 
result of expansion in the initial coupling g is denoted by oo. The arrows show points 
where A = Aq. 
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window /whole 

n = 2 

n = 1 

h = 0 

TO = —8 

0.993 

0.993 

0.961 

TO = —5 

0.940 

0.940 

0.908 


Table I. Ratio of the bound state eigenvalue of the small window hamiltonian with 
indices limited by rh and h, to the eigenvalue of the whole effective hamiltonian at 
A = 1 GeV calculated using expansion up to second power in the running coupling 
§(lGeV). 0.993 corresponds to the absolute accuracy of the bound state eigenvalue 
equal 12% and 0.908 to 19% (see the text). 


can project the narrow effective hamiltonian and reproduce the bound state 
eigenvalue by diagonalization of the projected matrix, is answered in Table I. 
The eigenvalue of the whole 7f(A = IGeV) is equal —0.8902 GeV instead of 
— 1 GeV. Window matrices with energy range order 1 GeV reproduce the same 
result with accuracy given in Table I. This is encouraging to pursue a similar 
approach to QCD. 

3. FOCK SPACE METHOD 

The model study shows that the perturbative similarity renormalization group 
allows a calculation of a small width effective hamiltonian, which can be pro¬ 
jected on a small space of states. The small hamiltonian can be solved exactly 
and the bound state eigenvalue of the full theory is obtained with 10% accuracy. 
The question is how to repeat these steps in quantum field theory. 

The method we propose 0 is based on the idea that one can unitarily trans¬ 
form the creation and annihilation operators, i.e. 

= Uxalul , (3.1) 

and the same for o’s. and aoo appear in the initial hamiltonian H. We call 
them “bare”. a\ and a\ appear in the effective hamiltonian TC\. They create 
and annihilate effective particles. In a way, U\ is analogous to the Melosh 
transformation in the case of quarks. However, we are building the transfor¬ 
mation using the similarity renormalization group idea, the transformation is 
fully dynamical and it can be applied to other particles than quarks, too. 

The effective hamiltonians satisfy the equation 

^nx = [nx,T^], (3.2) 

where T\ = u\dUx/d\. The unitary transformation generator T is constructed 
so that the effective hamiltonians have width A in the relative momentum 
transfer, 


'Hx = Fx[gx]. 


(3.3) 
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The operation F\ on the interaction terms Q\, inserts the similarity factors, f\. 
They are most easy to think about as form factors in the interaction vertices. 
The smaller is A the softer are the interactions and the effective particles get 
more dressed. 

Following the general idea of the similarity scheme one can hnd the 
equation satisfied by the vertex operators in the effective hamiltonians Q, i.e. 



/a02A, 1 “ f\)G2\ 


J Qix 


(3.4) 


Q\ = Gi\ + G 2 \, Gi\ is the a^o part of the hamiltonian and G 2 \ is the remaining 
part which changes momenta of the individual particles. The curly bracket with 
subscript Gi\ denotes the similarity energy denominator factor. 

An elementary example illustrates how it works in Yukawa theory which is 
defined by the following initial hamiltonian 


Hy = dx 




1 +-9-‘-^+TO^ 1 12 2 \, 

2id'^ “I" M 


2 - 7 

Fgip “h g 1pm n 


. a ;+=0 


The one particle energy is obtained in the form, 

Ql meson A — [^] 


k^2^,.2 


+ Ma t 

-aj.Ofe. 


k+ 


(3.5) 


(3.6) 


In second order perturbation theory in the coupling constant g, Eq. 
implies 




d/^(z^)8 (x-l)W 
dX A 42 _ ^2 


re{x,K). 


(3.4) 

(3.7) 


Here = (k^ + m‘^)/x{l — x) and re{x, k) denotes the regularization factor 
which is an analog of the number N in the matrix model. The similarity 
function P{z\) can be made as simple as, for example, 9[X^ + 3^^ — A4^). In 
this case, integration of Eq. (3.7) gives the following effective meson mass term 


z z I 

+ 4^ 


+ {y? - Qw?) log ^ 
^ 1 . 


+ -^i) + 0 ( 5 ^)- (3-8) 


Ai) denotes a finite term which has a limit when A ^ 00 . It equals 0 for 
A = Ai. fj-i is the effective meson mass in the hamiltonian ?-f(Ai). In the second 
order calculation, it is equal to the physical meson mass if Af < 4 to^ — 3/i^. 

The reason I show this example is that one can do similar calculations for 
other terms in the effective hamiltonians. @] Eor example, in second order 
perturbation theory, effective interactions between quarks are partly similar to 
the results obtained by Perry and his collaborators. @ @ 
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The questions how many orders of pertrubation theory are required in the 
calculation of the effective hamiltonian for constituent quarks and gluons in 
QCD and how large must be the subspace of the light-front Fock space to di¬ 
agonalize the effective hamiltonian of QCD, require much more work to answer 
than in the matrix model. 
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